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Abstract 



We present the design, the fabrication and the characterization of a 
tunable one-dimensional (ID) photonic crystal cavity (PCC) etched on two 
vertically-coupled GaAs nanobeams. A novel fabrication method which 
prevents their adhesion under capillary forces is introduced. We discuss a 
design to increase the flexibility of the structure and we demonstrate a large 
reversible and controllable electromechanical wavelength tuning (> 15 nm) 
of the cavity modes. 

1 Introduction 

The use of Nano Opto Electro Mechanical Systems (NOEMS) for the real-time 
spectral reconfiguration of a photonic crystal cavity (PCC) has recently drawn 
a lot of attention because of the large achievable tuning rates with small optical 
losses [1-3] and of the possibility of application in tunable filters [4] and cavity 
optomechanics [5]. Moreover, a tunable PCC coupled to quantum dots (QDs) 
facilitates the spectral alignment of the cavity to the emitter, opening up the 
opportunity to realize efficient, on-chip and scalable single photon sources [6] 
and to study cavity quantum electrodynamics phenomena [7]. Most of these 
systems are usually based on the electromechanical control over the evanescent 
coupling of two almost identical semiconductor cavities. In a previous work [H] 
we demonstrated a tunable InGaAsP double-membrane PCC whose resonant 
wavelengths could be electrostatically controlled over a 10 nm range. The use 
of vertically coupled cavities allows the simultaneous tuning of the cavity and 
the additional electrical control over active layers, located in one of the two 
membranes, e. g. by Stark tuning. This article describes the design, the 
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fabrication and the characterization of a tunable ID nanobeam cavity on GaAs 
with embedded InAs QDs. ID cavities have various advantages compared to 
2D cavities such as the smaller mode volumes, higher quality factors Q, higher 
compactness and ease of design [9j [10]. Electrostatically tunable ID PCCs 
have been already demonstrated using in-plane actuation of laterally-coupled 
nanobeams [21 13] while the vertically-coupled configuration and the electrostatic 
tuning of nanobeam cavities on GaAs has not been reported yet. We achieve 
electromechanical tuning over 15 nm with 15 V applied voltage. 

2 Design and theory of ID PC nanobeams 

The nanobeam PCC discussed in this work consists of a row of evenly spaced 
holes (lattice constant a) etched in a semiconductor beam of width w and rect- 
angular cross-section. By increasing one of the hole-to-hole spacing to 1.4a a 
cavity is formed. Such a simple design provides a single-mode cavity [11] and it 
has been chosen for its simplicity. Fig. []Ja) shows the in-plane TE-like mode 
profile (E y component) obtained by solving Maxwell equations in two dimen- 
sions numerically. In the configuration studied here, two nanobeams, one on 



Figure 1: (a) Geometry of the ID PCC on nanobeams and simulated in-plane mode 
profile of the E y component. The symmetric (left) and anti-symmetric (right) vertical 
profiles of the coupled system is also shown, (b) Sketch of the proposed diode struc- 
ture to realize short and tunable nanobeams: a 12 x 12/xm 2 doubly-clamped bridge 
with a 8/zm long nanobeam in the center. Only the holes and the side trenches are 
etched through both membranes. By operating the junction under reverse bias, the 
electrostatic force bends the upper slab and brings the nanobeams at a closer distance. 

top of the other, are brought at a very close distance (< 200 nm) to obtain op- 
tical coupling. This results in an energy splitting of the y-polarized mode into 
two modes having an anti-symmetric (at higher energy) and a symmetric (at 
lower energy) vertical profile (Fig. QJa)). By controlling the distance between 
the beams, the coupling strength can be modulated, resulting in a wavelength 
tuning of the cavity. To estimate the amount of coupling as a function of the 
geometrical parameters, a three-dimensional (3D) finite element method (FEM) 
is used [12] with a geometry adapted from [13] . By taking advantage of the sym- 
metries in the geometry, only one eighth of the double nanobeam structure is 
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simulated applying adequate boundary conditions on the symmetry planes. By 
enforcing a symmetric electric field (perfect magnetic conductor) or a symmetric 
magnetic field (perfect electric conductor) on the z-direction, the symmetric or 
the anti-symmetric modes can be calculated. A perfectly matched layer is used 
to simulate open boundaries. 3D simulations can be directly compared to pho- 
toluminescence (PL) experiments by integrating the radiated power (Poynting 
vector) from a dipole source [H] in the cavity to the surrounding air domain, 
assuming an ideal objective, capable of collecting light emitted in all directions. 
Fig. E(a) shows the calculated cavity resonant wavelengths as a function of the 
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Figure 2: (a) Results of a 3D FEM simulation of double nanobeams. The peak wave- 
length of the symmetric and the anti-symmetric modes are plotted against the inter- 
membrane distance. The parameters are: thickness t = 160 nm, width w = 420 nm, 
lattice spacing a = 370 nm and hole radius r = 96 nm. The refractive index dispersion 
in GaAs is taken into account. Ao is the wavelength in the uncoupled case (b) Calcu- 
lated Q factor, obtained from the solution of the lossy eigenvalue problem. Qo is the Q 
factor in the uncoupled case. 

nanobeam distance. A large shift (AA = 180 nm) from the uncoupled nanobeam 
case is achievable when the distance goes to zero. At the nominal distance cho- 
sen for the experimental realization of the device {zq = 200 nm), a maximum 
tuning before pull-in (i.e. at \z — zq\/ zq = 1/3) of AA = 23 nm and a tuning rate 
dX/dz = 0.2 nm/nm are predicted. The Q factor plotted in Fig. [2(b) is obtained 
solving an eigenvalue problem linearized around the mode wavelength. In the 
coupled region [z < 500 nm) the symmetric mode shifts to lower frequencies, 
away from the light cone, and consequently its Q increases. Conversely, the Q 
of the anti-symmetric mode drops almost to zero. The Q behavior, apparently 
in contrast with the double slab case p], is due to the Bragg mirror design 
(without tapers) which provides an abrupt termination of the cavity field and 
therefore a large amount of leaky components and radiation losses [9J [15] . As 
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the nanobeams are moved further apart, a modulation of the Q is observed even 
without mode coupling. This can be explained by interference effects which 
cause an amplification or a cancellation of the leaky modes [16] . 

To complete the NOEMS model, an electro-mechanical analysis of the device 
is performed. An important mechanical design parameter is the overall stiffness 
of the NOEMS represented by the Hooke's spring constant per unit area k 
(assuming a lumped electrostatic model made of metallic capacitor plates con- 
nected to springs). Taking the center of a doubly-clamped beam as the single 
degree of freedom of our model, the spring constant is given by k = 32Et 3 /L 4 
where E is the Young modulus {E = 85.9 GPa for GaAs), t and L are the 
nanobeam thickness (t = 160 nm) and length, respectively [17] . Assuming an 
uniform electric field, it is possible to derive the electrostatic force per unit area 
acting on each beam: P e iectrostatic = e o|p' where z is the distance between them 
and U is the total applied bias. Applying the equilibrium condition to the sym- 
metric system where two identical nanobeams are moving provides the equation 
^Peiectrostatic = k(zo — z), zq being the distance at rest. The displacement curve 
and pull-in voltages can be calculated solving: 

z 3 -z z 2 + e -^ = (1) 

A reasonable choice for k is such that pull-in occurs at U < 10 V, which is 
also a typical breakdown value for our p-i-n diodes. This restricts the stiffness 
of the nanobeams to k < 800 Pa/nm and consequently sets a lower bound to 
their length (L > 11/xm). As it will be discussed in section 3, it is possible to 
fabricate even longer structures. However a very large bending after undercut 
is observed (Fig. EJc)) leading to unpredictable coupling configurations and 
altering the mechanical properties of the nanostructure. To keep the nanobeams 
shorter and, at the same time, to lower the demand of actuation voltage in the 
devices, the nanobeams can be mounted on a larger and more flexible frame 
structure which allows to rigidly translate the upper beam to the fixed bottom 
one. The proposed device geometry is shown in Fig. H^b). Such a structure is 
not only more flexible but it also guarantees a more uniform application of the 
electric field than the simple nanobeam geometry. Due to the complexity of the 
shape, an effective stiffness k e fj, which can be calculated via FEM simulations 
or extracted from the measurements, is introduced. Since only the upper slab 
moves under the electrostatic pressure, the mechanical displacement corresponds 
to the distance between the nanobeams (z) and it is described by Eq. 1 replacing 
the stiffness by k = 2k e jj. To satisfy the conditions discussed above, an ideal 
value of k e ff should be lower than 400 Pa/nm. From FEM simulations, the 
structure stiffness ranges from 0.1 to 1 kPa/nm for the experimentally realized 
devices, depending on the actual geometry and dimensions. 
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3 Sample fabrication 



The sample is grown by molecular beam epitaxy on an undoped (100) GaAs 
substrate. A thick (1/xm) Alo.7Gao.3As sacrificial layer is initially deposited to 
separate the double-membrane structure from the substrate. Then, two GaAs 
membrane layers having the same thickness (160 nm) and an Alo.7Gao.3As inter- 
membrane layer (200 nm), as verified by scanning electron microscope (SEM) 
analysis, are grown. The upper membrane contains a layer of low density self- 
assembled InAs quantum dots, emitting around 1.3/im at room temperature 
[18j . To realize the electrostatic actuator, part of the membranes are doped to 
form a p-i-n junction. The top 50 nm of the lower membrane are p-doped and 
the bottom 50 nm of the upper membrane are n-doped (n = p = 3 • 10 18 cm -3 ). 
Since the QDs are situated above the junction, they are not affected by the 
electrostatic field. 

To fabricate the double nanobeams, a 400-nm-thick S13N4 hard mask is 
first deposited by plasma enhanced chemical vapor deposition (PECVD). The 
proximity-corrected nanobeam design, consisting of holes and side trenches, is 
patterned by a 30 keV electron beam lithography (RAITH-TWO 150) on a 360 
nm thick electron beam resist (ZEP 520A), aligned to the [011] or [0-11] direc- 
tions. After development, the nanobeams are first transferred to the underlying 
hard mask using pure CHF3 reactive ion etching (RIE) and then deeply etched 
(« 800 nm) through both GaAs membranes by inductively coupled plasma 
(ICP) (CI2/N2 chemistry at 200 °C). To remove the residual Si 3 N 4 mask and 
the sacrificial layer, a selective wet etching step in hydrofluoric acid is usually 
done at this stage. By doing so, however, the nanobeams will pin together dur- 
ing the drying step because of the very strong capillary forces developed and 
the low elastocapillary number (< 1) of the nanobeams in water |19j . In [8] we 
reported a method to avoid stiction without resorting to supercritical drying. 
It relies on the use of the S13N4 mask to stiffen the structure during the drying 
process. For nanobeams this method is not directly applicable because after 
ICP etching the mask supports only the top nanobeam leaving the other one 
free to collapse. Here we introduce a new technique which consists in the fabri- 
cation of nitride sidewalks around nanobeams before the undercut. The process 
steps are summarized in Fig. [3ja). 

The sample is cleaned in oxygen plasma and dipped into diluted HFiH^O 
(1:100) for 10 seconds to smoothen the bottom of the recesses without starting 
the AlGaAs undercut. Then a second deposition of a 600-nm-thick SisN4 con- 
formal layer is performed. The nitride is not deposited inside the holes, due to 
their small size (r < 100 nm). Using CHF3/O2 RIE at high power, the nitride is 
etched with a strongly anisotropic profile. Since the etch rate is much higher in 
the vertical than in the horizontal direction, a 300-nm-thick SisN4 supporting 
layer is left on the side of the nanobeams. By carefully optimizing the RIE 
times, the holes are opened again without damaging the GaAs. Fig. Gth) shows 
a SEM picture of the device cross-section after this step. Once the sidewalks are 
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fabricated and cleaned, no SisN4 is observed inside the holes and the sample is 
ready for the undercut. 



(a) ■ GaAs 




Figure 3: (a) The fabrication process used to realize freestanding nanobeams (seen 
in cross-section), (b) SEM picture of the cleaved hole cross-section after the sidcwall 
fabrication but before wet undercut. The sidewalls cover the nanobeam but the holes 
are opened again. No SisN4 is visible inside the hole, (c) 15 ^m long free-standing 
nanobeams. Stress relaxation induces a large bending in different directions which may 
also cause the nanobeams to touch and adhere. 

To remove the sacrificial layer selectively with respect to both GaAs and 
S13N4, a cold (1°C) HC1 36% solution is used [20] followed by a fast (3 seconds) 
dip in HF 5% to remove possible residues. The acid etches the AlGaAs through 
the top holes and the sides. Subsequently, the sample is rinsed in ultra-pure 
water and dried with nitrogen. The SisN4 sidewalls hold the nanobeams laterally 
and thus prevent their adhesion during drying. Finally, the nitride is removed 
isotropically in a low power CF4 plasma to minimize surface damage and to 
release the structures. The resulting free-standing nanobeams are shown in 
Fig. EIc). On long (> 10/rni) freestanding structures, the upper and the lower 
nanobeam relax and bend either upwards or downwards. This situation is not 
desirable, since it causes a non-reproducible coupling after fabrication and may 
also cause the nanobeams to stick together. For this reason shorter nanobeams 
are fabricated and placed on larger (hence more flexible) movable frame, as 
discussed above. 

To realize tunable structures, the process described above is realized on a 
previously prepared p-i-n junction with metal contacts. The diode's fabrication 
consists in two lithographic steps followed by wet etching (in citric acid/peroxide 
and HF 1%) to open vias to the p- and to the n-layer. Lift-off of Ti/Au 50/200 
nm is used to realize the contact pads. As discussed above, the larger structure 
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sketched in Fig. QJb), which serves as a supporting frame for the nanobeams, is 
defined during the p-via etching. All the movable parts are opened with holes 
to facilitate the undercut. The final device is shown in Fig. Eta). 

4 Measurement and results 

The device is tested with a micro-PL setup equipped with electrical probes. 
The QDs are pumped non-resonantly using a 785 nm diode laser from the top 
through a microscope objective (NA 0.4). The PL signal is collected from the 
objective, separated from the pump laser using a dichroic mirror, coupled into a 
fiber and analyzed with a spectrometer. The device is operated in reverse bias 
and for each voltage a spectrum is acquired. The measurements reported here 
have been obtained with the device of Fig. H^a). The actual geometry of the 
device has been measured by SEM. The results are summarized in Table HJ 

Table 1: PCC and nanobeam geometry as measured from SEM. When not explicitly 
specified, an uncertainty on the measurement of ±10nm is expected. 



Nanobeam thickness 


t 


150 nm 


Nanobeam width 


w 


420 nm 


Nanobeam length 


L 


(8 ± 0.2) fjcm 


PCC lattice constant 


a 


370 nm 


Hole radius 


r 


96 nm 



Fig. [4](b) shows the comparison between the PL and the FEM-simulated 
spectra using the parameters of Tableland different air gaps. A good agreement 
is observed for an initial gap zq = 175 nm and for the maximum displacement 
of 35 nm (z = 140 nm). The calculated initial gap is lower than the nominal 
thickness of the inter-membrane sacrificial layer because of the relaxation of the 
structure after undercut. As expected, the single-mode cavity shows a double 
peak due to the coupling of the nanobeams. The symmetric (s-)mode (initially 
at A = 1270 nm) and the antisymmetric (as-)mode (at A = 1190 nm) shift to 
longer and shorter wavelengths with increasing bias, respectively. The maximum 
shift is AA as = —15.4 nm for the as- mode and AA S = +13.6 nm for the s-mode. 
The tuning curve as a function of the applied reverse bias is plotted in Fig. 
Hl^c). The measured Q factor is Q s = 740 =L 40 for the symmetric mode and 
Q as = 450 ± 20 for the anti-symmetric mode. For the range of displacement 
considered here, the tuning of the Q, predicted by simulations, is not visible 
whereas the expected difference between the Q factors (Q s > Q as ) is observed. 
The low Q compared to simulations is attributed to the poor selectivity of CF4 
plasma towards GaAs which damages the holes and reduces the thickness of the 
upper nanobeam. The latter is also causing the small difference of tuning rates 
and PL intensity between the s and the as mode. In the PL spectra, two quasi- 
degenerate band-edge modes (matching the simulated data) are also visible at 
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Figure 4: (a) SEM image of the final device used for the tuning experiments, (b) PL 
spectra at and 15 V DC bias (maximum tuning) compared to FEM simulations, (c) 
Tuning curve of the observed ID PCC peaks, showing the as-mode, the s-mode and a 
double antisymmetric dielectric band-edge mode, (d) Calculated displacement of the 
upper membrane as a function of the voltage. The curve has been fit with a lumped 
electrostatic model to extract the equivalent device stiffness /c e //. 
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A = 1326 nm at V. These modes arise from slow-light PL enhancement |21] 
at the dielectric band edge and they are de-localized over the nanobeam length. 
Due to the double layer structure, they also split into s and as profiles. A higher 
tuning range compared to the cavity (AXbe = —21.2 nm) is obtained. This can 
be explained considering that the electromagnetic field of the dielectric BE has 
an in-plane distribution which is mostly located in the high-index area of the 
GaAs nanobeam, therefore modulation of the effective index in these areas is 
expected to affect these modes more than the cavity mode. This has been also 
verified by 3D band calculations of the double photonic crystal nanobeam and 
the resulting shift of the anti-symmetric dielectric band has been plotted in Fig. 
S(b) (dashed line). 

From the value of the coupling it is possible to estimate via FEM simulations 
the distance between the nanobeams and to derive the displacement as a func- 
tion of the applied bias (Fig. H^c) squares). The curve has been fitted using the 
lumped model of Eq. 1 (in the case of fixed bottom membrane) to estimate the 
effective stiffness k e ff of the entire device. A value of k e jf = (1.1 ± 0.1) kPa/nm 
is obtained. This value is still higher than the design parameters discussed in 
Section [2] and more flexible designs are needed to achieve even higher tuning 
ranges. The displacement curve also shows a large deviation from theory when 
the bias voltage is > 11 V. The saturation is due to the breakdown of the p-i-n 
junction, causing the currents in the intrinsic region to increase and limiting the 
voltage across the air gap. 

5 Conclusions 

A GaAs NOEMS device for the wavelength tuning of a ID PCC on nanobeams 
has been demonstrated. By designing a flexible structure, a large reversible shift 
of the anti-symmetric mode (AA as = 15.4 nm) has been observed with an ap- 
plied bias of U = 15 V. A good agreement with FEM simulations has also been 
obtained. Higher Q cavities may be obtained applying more sophisticated PCC 
designs and improving the fabrication process. The ID nature of nanobeams 
opens up several opportunities for the design of on-chip tunable filters. It also 
enables the realization of a mixed in-plane and out-of-plane actuator, combining 
up to four nanobeams, to increase the tuning even further. Moreover, the pos- 
sibility to combine a tuning control with an active region, makes such a device 
highly attractive for the realization of tunable lasers, single photon sources and 
quantum photonic integrated circuits. 
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